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Abstract 

We analyze several continuum models of polymers: worm-like chains, ribbons and Fourier knots. 
We show that the torsion of worm-like chains diverges and conclude that such chains can not be 
described by the Frenet-Serret (FS) equation of space curves. While the same holds for ribbons 
as well, their rate of twist is finite and, therefore, they can be described by the generalized FS 
equation of stripes. Finally, Fourier knots have finite curvature and torsion and, therefore, are 
sufficiently smooth to be described by the FS equation of space curves. 
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I. INTRODUCTION 



Recent progress in the ability to manipulate single biomolecules such as double stranded 
DNA and protein filaments , prompted the development of continuum models of complex 
polymers capable of describing bending fluctuations and finite extensibility under tension 
(extensions of the worm-like chain model ^) as well as models that can describe twist rigid- 
ity, spontaneous twist and chiral response to torque (the ribbon model [3]). These models 
are defined by their elastic energy functions which are then used to generate the equilibrium 
ensemble of polymer conformations, based on the conventional Gibbs distribution approach 
(i.e., weighting the conformations by an appropriate Boltzmann factor). While this approach 
proved to be quite successful for open polymers , it could not be apphed to generate the 
conformations of polymer loops and to study properties of circular double stranded DNA 
such as supercoiling and formation of knots. In order to cope with the latter problem, we 
developed a purely mathematical procedure of generating closed curves, based on the expan- 
sion of the components of the polymer conformation vector f{t) {t is some parametrization 
of the contour of the loop) in finite Fourier series and taking the corresponding Fourier 



coefficients from some random distribution 



6|. We found that this distribution of Fourier 



knots could be fine tuned to mimic some of the large scale properties of closed Gaussian 
loops and small scale properties of worm-like chain, this could not be achieved using a single 
persistence length. 

The present deals attempts to establish a common framework for the discussion of the 
above physical and mathematical models of polymers, by classifying them according to 
their smoothness. A continuous curve f{s) {s is the contour parameter which measures the 
distance along the contour) is defined as an m-smooth curve if its mth derivative d'^f/ ds^ 
is a continuous function of s. In section [Til we introduce the fundamental equations of 
differential geometry of space curves and of stripes and show that these equations describe 
space curves that are at least 2-smooth and stripes that are at least 1-smooth. This is 
equivalent to saying that while such space curves must have finite torsion, stripes must 
have finite rate of twist but their torsion is free to take any value along their contour. In 
section UTTl we show that worm-like chains belong to the class of freely rotating models, with 
uniformly distributed dihedral angles, divergent torsion and a normal whose direction jumps 
discontinuously as one moves along the contour of the chain and, therefore, such objects can 
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not be described by the Frenet-Serret (FS) equation. We also show that because the energy 
of a ribbon depends on its twist, the typical conformations of ribbons have finite rate of 
twist but their torsion diverges at many points along the contour. In section HVl we compare 
the ensembles of 1-smooth (ribbons) and oo-smooth (Fourier knots) curves and show that 
typical realizations of the latter (but not the former) ensemble, have finite torsion and a 
smoothly varying normal and can be described by the FS equation. We also calculate the 
distributions of spatial distances between two points on the contour of the curve in the above 
ensembles and find that these distributions differ significantly only for distances of the order 
of persistence length. Finally, in section |V] we discuss our results and conclude that unlike 
worm-like chains and ribbons which possess no torsional rigidity, the ensemble of curves 
generated by the Fourier knot algorithm can be characterized by finite torsional persistence 
length. 



II. DIFFERENTIAL GEOMETRY OF CURVES AND STRIPES 



It is often convenient to represent a space curve r(s) defined in a space-fixed coordinate 
frame by intrinsic coordinates, as follows. At every point along the curve (0 < s < L) one 
constructs a set of three orthogonal unit vectors known as the Frenet frame: the tangent 
i defined as i(s) = df/ds, the normal n(s) which points in the direction of di/ds and the 
binormal h{s) = t{s) xn{s). The rotation of the Frenet frame as one moves along the contour 
of the curve is described by the Frenet-Serret (FS) equation [tI: 
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where k is the curvature and r is the torsion (in general, both are functions of s). The 
condition of validity of the above equation is that nds, rds — > as (is — * everywhere along 
the curve. It is straightforward to show that 
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and we conclude that since the condition of validity of the FS equation is that the torsion 
is finite everywhere along the curve, the curve should be at least 2-smooth. 
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Given the curvature and the torsion at each point along the curve, one can solve Eq. ([T]), 
calculate the tangent t{s) and integrate it to construct the parametric representation of the 
space curve, r(s). The above construction is unique in the sense that any pair of functions 
K,{s) and t{s) defined on the interval < s < L, can be uniquely mapped to a space curve 
r(s) of length L. The simplest examples are (a) k, = const., t = which yields a planar 
circle of radius and (b) k = const., r = const, which corresponds to a helix. 

We now turn to consider stripes of length L, width W (such that L ^ W ) and thickness 
D — > 0. Unlike a space curve which is uniquely defined by the tangent vector i{s) (the normal 
and the binormal are auxiliary constructs, needed only to calculate the tangent, given the 
curvature and the torsion), a stripe is a slice of a plane with which one can associate two 
orthogonal unit vectors ti (in plane) and ^2 (normal to the plane). The spatial configuration 
of the stripe is thus defined by the local orientation (at each point s on the centerline) of 
the orthogonal triad known as the Darboux frame which specifies the directions of the two 
axes ii and £2 and that of the tangent to the centerline that runs along the long axis of the 
stripe, i = £3. While a space curve is completely defined by the two functions k{s) and r(s), 
a stripe is represented by three generalized curvatures ujk{s) (k = 1,2, 3) that determine the 
unit vectors {tj} via the generalized Frenet-Serret equation 
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This can be written compactly as tj = Ai^ik where Ai^ = —tijkUJj {tijk is the Levi-Civita 
tensor). Inspection of Eq. ^ shows that uji{s)ds is the infinitesimal angle of rotation 
about the direction ti{s) and the condition of validity of the generalized FS equation is that 
this angle vanishes in the limit ds 0. However, unlike the torsion r which is completely 
determined by the space curve f{s) and can be expressed in terms of its first three derivatives, 
the rate of twist is the local rate (per unit length) of rotation about the tangent to this 
curve and, as such, it depends only on the orientation of ti and t2, and can not be expressed 
in terms of the centerline r(s) and its derivatives! We conclude that in order for Eq. ([3D to 
hold, the centerline of the stripe should be represented by a 1-smooth curve. 

Since both couples of unit vectors jn, foj and {^1,^2} lie in the plane perpendicular to 
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the local tangent, the triads |t, n, fej and {^3,^1,4} are connected by 
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where the matrix 



1 

S = COS a sin a (5) 
— sin a cos a j 

generates a rotation by an angle a about the i = axis. Assuming that the centerline of 
the stripe can be represented by the FS equation (i.e., that it is an 2-smooth curve), one 
can express the generalized FS equation in terms of the curvature, torsion and the angle a 
between the directions of the binormal h and the ti axis: 
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Comparing ([6]) with ([3]) yields: 
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Conversely, 



da/ds 
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(10) 



r = 0^3 — da/ds 

However, since the centerline of the stripe is required to be only 1 -smooth, these relations 
do not hold in general! 



III. POLYMER MODELS: WORM-LIKE CHAINS AND RIBBONS 



Physical models of polymers are based on the choice of a geometrical model and an energy 
functional. The simplest and the most prevalent model is that of a continuous Gaussian 
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random walk, with which one can associate a free energy that describes the entropic cost 
of stretching the polymer chain, Egrw = \cikBT j^^ dt {dr/ dtf . Here a is the "monomer" 
(cutoff) length, is the Boltzmann constant and T is the temperature j^. Notice that this 
free energy is expressed only in terms of the first derivative of the trajectory, dr/dt and, 
therefore, space curves that describe polymer conformations in the continuous Gaussian 
random walk model have to be only 0-smooth (only the curve itself and not its derivatives, 
has to be continuous everywhere). 

The worm-like chain model of polymers combines bending elasticity and inextensibility (the 
latter condition can be expressed as \dr/ds\ = 1) and, assuming that the stress-free state 
corresponds to a straight line, the energy can be written as |lol | 
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(11) 



Since the energy depends only on the curvature k and does not depend on the torsion r, 
the corresponding space curve has only 1 -smooth. Note that curves described by the FS 
equations have to be at least 2-smooth and, therefore, worm-like chains are not sufficiently 
smooth to be described by the fundamental equations of differential geometry of space curves! 
In order to get physical intuition about the origin of the problem, lets consider a discretized 
model of a continuous curve in which the polymer is made up of connected straight segments 
of length As each, such that the direction of the segment at point s is given by the tangent 
to the original chain at this point, i{s) (the continuum limit is recovered as As 0). The 
angle between neighboring segments is denoted as A^(s) and, in order to describe the non- 
planar character of a general space curve, one has to introduce the dihedral angle A<^(s) 
between the two successive planes {t(s), t(s + As)} and [i{s + As), t(s + 2As)} determined 
by three successive segments at points s, s + As and s + 2As. Since Ay9(s) is also the angle 
between neighboring binormals h{s) and 6(s + As), the Frenet frames at points s and s + As 
are related by a simple rotation 
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where the rotation matrix B [s) is given by 
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cos A9 sin A9 

- cos A(/9 sin cos At/) cos A^^ sinA^j 
sinAy9sinA^^ — sin Ay9 cos A^^ cosAy9 

Notice that no assumption is made so far about the magnitude of the angles A6 and AyD. 
Defining Ai{s) = i{s + As) - i{s), A h{s) = n{s + As) - n{s) and Ab{s) = b{s + As) - b{s) 
and subtracting the vector (i, n, b) from both sides of Eq. (fT3l) . yields 
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with / the unit matrix. Notice that unlike the FS equation which is valid only for infinites- 
imal rotations of the Frenet frame, Eq. f[T^ describes finite rotations; the FS equation can 
be derived from it by dividing both sides of the equation by As and taking the limit As 
( lim At/ As = dt/ds, etc.). In order for the right hand side of Eq. (IT^ to remain finite 
in this limit, all the elements of B (s) — I have to vanish .Since this is equivalent to the 
condition A9, Aip <^ 1 , one can expand the cosine and sine functions in Eq. ([13]) and, upon 
substituting the result into Eq. (IHj), one recovers the FS equation with k, = hm AO / As 
and r = lim Aud/As. 

As^O 

Returning to the worm-like chain model, we notice that while the bending energy ensures 
that the curvature is finite and the angle A6 is always small, there is no corresponding 
physical restriction on the magnitude of A(p and we conclude that the worm-like chain 
corresponds to the class of freely rotating chain models in which the angle Aip can attain 
any value in the interval [— tt, tt] . For such models, the expansion of B (s) — I in terms 
of Aip breaks down and the corresponding (1-smooth) curves can not be described by the 
FS equation. Notice that if one keeps the definition r = hm Aip/ As, the torsion can 
diverge at any point along the curve, generating abrupt jumps of the normal (see Fig. [1] 
b). Nevertheless, since the shape of a curve is completely characterized by (and only by) 
the tangent to it and since the latter changes continuously even for 1-smooth curves, such a 
curve appears (to the eye) to be just as smooth as an cxD-smooth one. 

Let us now consider the ribbon model of polymers which was designed to take into account 
rigidity with respect to twist and spontaneous twist of complex polymers such as double 



stranded DNA. In general, the ribbon has an asymmetric cross section with symmetry axes 
ti and t2 and a centerhne described by the tangent, t^, such that the triad of unit vectors 
can be associated with the Darboux frame famihar from the differential geometry of 
stripes. In the framework of the linear theory of elasticity of slender rods [hI]), the energy of 
a particular configuration of a ribbon is a quadratic functional of the deviations of its three 
curvatures {ui{s)} from their equilibrium values in the stress-free state, {uQi{s)} [3|]: 



(15) 



^R = llo [^1 (^1 " ^oi) + h iuj2 - UJ02) 

Here bi and 62 are the bending rigidities associated with the two principal symmetry axes 
of the cross section and 63 is the twist rigidity (the persistence lengths {aj} are obtained by 
dividing the corresponding rigidities by ksT). For ribbons with a symmetric cross-section 
(61 = 62) and without spontaneous curvature and twist (cuoi = 0), the above expression can 
be simplified 
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where = ujI+uj"^- Since the above energy functional depends only on the second derivatives 
of the curve (curvature) and on the twist, and does not explicitly depend on the third deriva- 
tives (torsion), the torsion of the curve is not controlled by the elastic energy and, therefore, 
the centerline of the ribbon can not be described by the FS equation. To demonstrate it, 
we rewrite Eq.(l3]) in a discretized form, using (fT3l) and the relation ([5]): 
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Comparing (fT7|) to ([3]) gives, in the limit As 0: 



ujj, = Aa/As + Aip/As da/ds + t 

LUi = (A0/As) sin(a;) — > nsma (18) 
LJ2 = {A9 / As) cos{a) ftcosa 

The fact that the elastic energy of a ribbon introduces a penalty for bending and twist 
deformations, ensures that only configurations with cjjAs — » (i = 1, 2, 3) contribute in the 
continuum limit As — > and, therefore, the conformations of a ribbon can be described by 
the generalized FS equation ([3]). Inspection of Eq. (fT8l) shows that the condition that the 
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twist accumulated over a contour distance As, is small, can be expressed as Aa + A(y9 << 1. 
Since there is nothing that restricts the magnitudes of Aa and of Ay9 separately, they can 
be arbitrarily large provided that the condition Aa ^ —Aip is satisfied and, therefore, the 
torsion associated with the centerline of the ribbon can be arbitrarily large (recall that 
the torsion is defined as the limit of A(f/As as As 0). Indeed, inspection of a typical 
conformation of a ribbon (taken from the ensemble of conformations generated using the 

n 

algorithm described in ref. p) shows that even though the rate of rotation of the physi- 
cal axis ti is everywhere finite (Fig. [1^), the rate of rotation of the normal fi is not (Fig. [T)d)! 



IV. ENSEMBLES OF 1-SMOOTH AND oo-SMOOTH CURVES 

Increasing the degree of smoothness from m to m + 1 acts as a constraint that prohibits 
certain configurations of a curve, and it is interesting to compare the properties of curves 
with different degrees of smoothness. Such a comparison is meaningful only in a statistical 
sense and in the following we will consider some physically relevant statistical properties of 
worm-like chains for which analytical results are available, with those of computer generated 
ensembles of ribbons (1-smooth) and Fourier knots (oo-smooth). In order to generate the 
ensemble of centerlines of ribbons (1 -smooth curves), we use the so called "Frenet algorithm" 
described in detail in ref. . In view of the discussion in the preceding section, such curves 
can not be described by the FS equations and, in order to avoid possible misinterpretation, 
we will refer to it as the ribbon algorithm. 

Let us compare the distribution of the angle A(f between neighboring binormals (or rather 
of cos Ay9 = 6(s) ■ b{s + As)) in the discretized version of the worm-like chain model, with 
that obtained by generating the ensemble of ribbon conformations, computing the centerline 
of each conformation and extracting the distribution of cos{A(p). In the worm-like chain 
model Aif is distributed uniformly in the interval [— vr, it] and, therefore, the probability 
distribution of cos(A(/)) is given by 

PwLC (cos Aip) = — ^ =d (cos Av9) (19) 
tta/I — COS"' A(p 

Using the ribbon algorithm to obtain the ensemble of conformations of a ribbon with a 
symmetric cross section and without spontaneous curvature, we generate the corresponding 
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distribution Pr [cos(A(y9)]. Up to numerical accuracy we find that the above distribution 
coincides with the simple worm-like chain expression, Eq. [19] and does not depend on the 
bending or twist rigidity (see Fig. [2]). This concurs with our expectation that, just like 
worm-like chains, centerlines of ribbons can be described by freely rotating type models. 
In order to generate oo-smooth curves we use the Fourier knot algorithm which was originally 
developed with the goal of investigating knots (i.e., closed curves) [6;]. Unlike methods 
Dased on modeling the knot as a 0-smooth curve made of discrete, freely jointed segments 



I2I . [13I . [l3|). this algorithm generates infinitely smooth knots, such that the derivatives 
\d^f/dt'^\ are finite for all m. The Fourier knot algorithm is based on the fact that for any 
closed curve parameterized by some arbitrary parameter t, the projections of the position 
vector f{t) on the Cartesian coordinate axes are periodic functions rj(t) = rj(t + T) with 
period T and can be expressed as finite Fourier sums: 



.w = E 



A>os(^) + 5;sm(^^ 



(20) 



n=l 

Different realizations of closed curves can be generated by choosing coefficients A^, i?^ from 
some statistical distribution. When the coefficients are given by An~^ exp — n/no, where A 
are random numbers in the interval [—1,1] and uq is an effective cutoff [riQ <^ ^max), the 
long wavelength properties of the ensemble generated by the Fourier knot algorithm, are 
in good agreement with those obtained from the worm-like chain model. These proper- 
ties include second moments such as the mean square distance between two points on the 
contour ([r(si) — r(s2)]^) and the tangent auto-correlation function ^ [t{si) ■ t{s2)] ) where 
\si — S2I >> /, with the persistence length I determined by the cutoff as / = 0.58L/nQ. How- 
ever, even though the tangent auto-correlation function decays exponentially with |si — S2I 
on length scales comparable to I (just like in the worm-like chain model), the correspond- 
ing decay length la = 0A35L/nQ is smaller than the persistence length obtained from the 
long-wavelength properties of Fourier knots. In ref. 6[] we suggested that the ensemble of 
configurations generated by the Fourier knot algorithm is equivalent to a physical ensemble 
of polymers which possess both bending and twist rigidity and, while the short range prop- 
erties of the tangent-tangent correlation function are determined by the bending persistence 
length only, both bending and twist persistence length control its long distance behavior. 
In any case, the fact that the ensemble of Fourier knots can not be characterized by a single 
persistence length suggests that the statistical properties of this ensemble differ from those 
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of worm-like chains and that it is important to investigate not only the second moments but 
the entire distributions. 

The first property we examine is the probability distribution Ppx (cos A;/)). As can be seen 
in Fig. [21 the distribution has a peak at cos At/) = 1, i.e., at Ayj = 0. Since A(y9 = rAs, 
we conclude that the ensemble of curves generated by the Fourier knot algorithm is char- 
acterized by finite torsion and a normal whose direction varies smoothly along the contour 
of each curve and, therefore, such curves can be described by the FS equation. We would 
like to stress that even though the torsion is described by the first three derivatives of r all 
of which are finite for Fourier knots, the observation that the ensemble of Fourier knots is 
dominated by curves with finite torsion is non-trivial since the expression for the torsion di- 
verges at points along the contour where the curvature vanishes (see Eq. H]). Notice that for 
ribbons with no spontaneous curvature and twist, the partition function can be written as 
the product of bending and twist parts Z = ZhendZtwist, with Z^end is given by the functional 
integral (assuming a symmetric ribbon of bending persistence length /) 



Since ujf + 002 = k^, the measure D {uji} D {0^2} can be written as the product of a "radial" 
contribution D{k}k, and an angular one. We therefore conclude that the probability of 



finite everywhere, as observed. Strictly speaking the above argument was derived for open 
ribbons and not to closed curves, but since it involves only the measure and not the form 
of the energy function, it applies to Fourier knots as well (see Fig. [3] where the measured 
distribution is plotted for Fourier knots with uq = 150). 

We now turn to compare the statistical properties of 00-smooth curves generated by the 
Fourier knot algorithm and the 1-smooth centerlines of ribbons generated by the ribbon algo- 
rithm. Consider the probability distribution P{R\sq) of the distance R{so) = \r{si) — r{s2)\ 
between points Si and S2 {sq = \si — S2I) along the contour (the second moment of this 
distribution for Fourier knots was calculated in ref. Q]). The difficulty in comparing the 
two ensembles is that while the ribbon algorithm generates open curves, the Fourier knot 
algorithm yields closed loops. In order to compare the latter with the former, we make 
use of the fact that, as long as we consider contour distances (sq) much shorter than the 





points with k 



vanishes linearly with k, and since r oc l/zt' 



the torsion should be 
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total length of the loop (L), Pfk{R\so) approaches the probability distribution for an open, 
infinitely smooth curve. In Fig. H] we plot Pfi{R\so) (ribbon) and Pfk{R\so) (knot). As 
expected, in the long wavelength limit (sq ^ I) the two distributions approach the Gaussian 
random walk result, Pgrw{R\so) -R^ exp[— 3-R^/(4so/)] (see blue triangles). Since on very 
short length scales {sq << I) all distribution functions approach the trivial limit 6{R — Sq), 
the two distributions can only differ on intermediate length scales (sq ~ I). This is indeed 
confirmed by our simulation results. Fig. HI Notice that in this regime the maximum of 
Pfk{R\so) is shifted to higher values of R than that of Pfi{R\so) indicating that typical 
conformations of 1-smooth curves are more compact than those of oo-smooth ones. The 
origin of the difference can be traced back to the fact that the characteristic magnitude of 
the torsion of a 1-smooth curve is much larger than that of an cxo-smooth one and, therefore, 
on length scales comparable to the persistence length, the latter curves are confined to a 
plane while the former have a three dimensional character. 



V. DISCUSSION 



We have demonstrated that the standard continuum models of polymers including con- 
tinuous Brownian random walks, worm-like chains and ribbons, generate space curves that 
are not sufficiently smooth to be described by the fundamental FS equation of differential 
geometry. Examination of the corresponding statistical ensembles shows that the dihedral 
angle Ay? between two successive binormals along the chain contour is uniformly distributed 
in the interval [— tt, tt] and we conclude that both worm-like chains and centerlines of ribbons 
belong to the class of freely rotating models, with divergent torsion and discontinuous jumps 
of the normal to the curve. However, unlike worm-like chains, ribbons have twist rigidity 
which means that the rate of twist of the physical axes of the cross section remains finite 
everywhere along the contour of the ribbon and guarantees that the triad of unit vectors 
associated with the ribbon obeys the generalized FS equation familiar from the differential 
geometry of stripes. We compared some statistical properties of ensembles of 1-smooth and 
oo-smooth curves generated by the ribbon and the Fourier knot algorithms, respectively. We 
showed that in the latter case the dihedral angle is peaked about Aip and, therefore, 
typical configurations of Fourier knots have finite torsion everywhere and can be described 
by the FS equation. We also compared the distribution functions of the spatial distance 
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between two points along the contour of a ribbon and of a Fourier knot. As expected, both 
distribution functions approach the hmiting Gaussian distribution for length scales much 
larger than the persistence length, but are quite different on length scales comparable to the 
persistence length. 

Finally we would like to stress that while the physical ensembles of conformations of worm- 
like chains and ribbons are generated using the standard methods of statistical physics (each 
conformation is weighted with an appropriate Boltzmann factor, exp{—E/ki)T)) , the ensem- 
ble generated by the Fourier knot algorithm is a purely mathematical construction and there 
is no elastic energy associated with different conformations of Fourier knots. Nevertheless, 



the observation of two persistence lengths reported in ref. {g] and the present finding that 
Fourier knots have finite torsion, suggest that the statistical properties of this mathemat- 
ical ensemble (notice that persistence lengths can be measured directly from the ensemble 



of conformations of the space curves, just as is done in AFM experiments [15|]) are quite 
similar to those of a physical ensemble of conformations of polymers with both bending and 
torsional rigidity. The detailed exploration of this analogy is the subject of future work. 
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FIG. 1: A typical conformation of a ribbon: the centerline is shown by the dark sohd curve and the 
gray hnes orthogonal to it show the direction of (a) one of the symmetry axes of the cross section 
ti and (b) the normal h 
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FIG. 2: (color online) The distribution function of cos{<p) of a worm- like chain (black solid line) 
and of a symmetric ribbon with ai = 02 = 0.1, 03 = 0.01 (red); ai = 02 = 0.1, = 1000 (green); 
ai = a2 = 0.01, 03 = 1000 (blue). The distribution of a Fourier knot with no = 150 is shown by 
the orange line. 
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FIG. 3: The distribution of curvature of Fourier knots. The sohd hne is a hnear fit to the distri- 
bution at K — 5- 0. 




FIG. 4: (color onHne) P{R\sq) vs R/sq is plotted for Fourier knots (full symbols) and symmetric 
ribbons (empty symbols), of total length L = 47r and contour distance sq = L/100. Taking the 
persistence length of the ribbon as I = 0.58*L/no we plot the distribution for: no = 50, I/sq = 1.16 
(black squares); no = 150, 1/sq = 0.386 (red circles); no = 600, I/sq = 0.096 (blue triangles). The 
Gaussian distribution with I/sq = 0.096 is shown by the solid cyan line. 
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